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Abstract. Critical exponents are computed for a variety of twist- 2 composite operators, which 
occur in polarized and unpolarized deep inelastic scattering, at leading order in the 1/Nt expan- 
sion. The resulting d-dimensional expressions, which depend on the moment of the operator, are 
in agreement with recent explicit two and three loop perturbative calculations. An interesting 
aspect of the critical point approach which is used, is that the anomalous dimensions of the 
flavour singlet eigenoperators, which diagonalize the perturbative mixing matrix, are computed 
directly. We also elucidate the treatment of 7 s at the fixed point which is important in simpli- 
fying the calculation for polarized operators. Finally, the anomalous dimension of the singlet 
axial current is determined at 0(1/Nf) by considering the renormalization of the anomaly in 
operator form. 
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1 Introduction. 



Our understanding of the structure of nucleons is derived primarily from experiments where 
they are bombarded by other nucleons or electrons at high energies. These deeply inelastic 
processes are, in general, well understood in most instances. Current activity, however, centres 
on examining polarized reactions due, for example, to the discrepancy observed in results for 
the spin of the proton and theoretical predictions by the EMC collaboration, Q. Consequently 
in order to make accurate statements about the data current theoretical interest has focussed 
on carrying out higher order perturbative calculations in the underlying field theory, quantum 
chromodynamics, (QCD). As this theory is asymptotically free at high energies, M, the coupling 
constant is sufficiently small so that perturbative calculations give a good description of the 
deep inelastic phenomenology. Indeed unpolarized scattering is well understood with one and 
two loop results available for the anomalous dimensions of the twist- 2 flavour non-singlet and 
singlet operators which arise in the operator product expansion, [3-5]. Moreover, the Dokshitzer, 
Gribov, Lipatov, Altarelli, Parisi, (DGLAP), splitting functions, [||, which are a measure of 
the probability that a constituent parton fragments into other partons are known to the same 
accuracy, ||, ||. The moments of the scattering amplitudes have also been studied. More recently 
the three loop structure has been obtained exactly, for low moments, through huge impressive 
analytic computations, both for non-singlet and singlet cases, 0. 

The situation for polarized scattering is less well established. The one loop anomalous 
dimensions for the corresponding twist- 2 (and 3) operators were computed by Ahmed and Ross 
in [||. However, only recently has the two loop structure been determined for the flavour 
singlet operators, B. (The non-singlet polarized dimensions are equivalent to the non-singlet 
unpolarized case.) This was checked by Vogelsang in [10, 11] by calculating the splitting functions 



themselves and then comparing with || by taking the inverse Mellin transform. In this the 
moment n of the operator has the conjugate variable x which is the momentum fraction of the 
parton in the nucleon. These results have been important for the next to leading order evolution 
of the structure functions to low x and Q 2 regions, p2| . (For review articles see fl3fl .) Crucial 
in this exercise is the dependence of the results on the moment of the operators. 

To go beyond this two loop picture would require a great amount of computation based, for 
example, on the unpolarized results of J?]]. One way of improving our knowledge would be to 
compute the appropriate quantities using another approximation. For example, the large Nf 
expansion, where Nf is the number of quark flavours, has been used to obtain the leading order 
coefficients of the anomalous dimensions of the twist-2 unpolarized non-singlet operators to all 
orders in the perturbative coupling constant, [jlj]. The resulting analytic function of n provided 
a useful check on the exact low moment non-singlet three loop calculation of 0. Briefly the 
method involved studying the scaling behaviour of the operator at the non-trivial fixed point in 
(i-dimensional QCD. With that success and the current interest in polarized physics required for 
exploring new x regimes it is appropriate to apply the large Nf analysis to study the dimensions 
of the underlying (twist-2) operators. Aside from providing n-dependent results and getting a 
flavour of the structure beyond two loops, it will give at least another partial check on the recent 
results of [9-11]. In particular we will determine the critical exponents at 0(1/Nf) which encode 
the all orders coefficients of the twist-2 polarized singlet operators. As a prelude we need to 
study the unpolarized singlet case which will extend the result of [Q. 

The paper is organised as follows. The basic formalism and notation is introduced in section 
2. We review previous work in section 3, including the technical details of the computation 
of singlet unpolarized operators anomalous dimensions. As the treatment of four dimensional 
objects, such as 7 5 , is needed we review previous 1/Nt work involving this in section 4. The 
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remaining sections 5 and 6 are devoted to the application of the results of earlier sections to 
polarized operators. In particular the latter section centers on the treatment of the singlet axial 
current which is not conserved due to the chiral anomaly. Future work and perspectives are 
discussed in section 7. An appendix gives details of the relation of the 1/Nf exponent results to 
the DGLAP splitting functions. 

2 Background. 

To begin with we review several of the more field theoretic aspects of the formalism including the 
role of the critical renormalization group. First, we recall that critical exponents are fundamen- 
tal quantities. In experiments and condensed matter problems dealing with phase transitions 
they completely characterize the physics. In order to describe such phenomena one determines 
estimates of the exponents by calculating, for example, in the underlying quantum field theory 
describing the transition. In practice this means carrying out a perturbative calculation of the 
renormalization constants of the theory, in some renormalization scheme. This information is 
then encoded in the corresponding renormalization group functions such as the anomalous di- 
mensions of the field or the mass. The critical coupling, g c , is subsequently determined from the 
/^-function of the theory. It is defined to be a non-trivial zero of /3(g). The appropriate critical 
exponents are then found by evaluating the anomalous dimensions at the critical coupling. In 
practice provided enough terms of the series have been computed, relatively accurate numerical 
estimates can be obtained. (Useful background material can be found in, for example, fll5| .) 

For the present problem we will examine a fixed point in QCD in d-dimensions and obtain the 
critical exponents as a function of d and A^ which characterize the transition. As indicated these 
correspond to the renormalization group equation, (RGE), functions evaluated at g c . Therefore 
provided the location of g c is known in some approximation like 1/Nf one can decode the 
information contained in the exponent and determine the coefficients of the RGE functions for 
non-critical values of the coupling. In this large Nf method it turns out that the structure of 
[3(g) is such that for leading order calculations in 1/Nf only knowledge of the one loop coefficient 
is required, [lq| . We illustrate these remarks with a general example, which will set notation 
for later sections. As we are interested in the coefficients in the series of an RGE function we 
denote such a function by 7(5) and define its expansion to be, with explicit Nf dependence, 

7(9) = a ig + (a 2 N f + b 1 )g 2 + (a 3 N f 2 + b 2 N f + Cl )g 3 + 0(g 4 ) (2.1) 

where the obvious definition for the 0(g 4 ) term is understood. The coefficients {aj, 6j, Cj, . . .} 
can of course be functions of other parameters such as colour group Casimirs or n. To evaluate 
(2.1) at a fixed point we take the general structure of the /3-function to be, in (i-dimensions, 

0(g) = (d - A)g + AN f g 2 + (BN f + C)g 3 + 0(g A ) (2.2) 

Setting d = 4 — 2e, then there is a non-trivial fixed point, g c , at 

* = + (wj) < 2 ' 3 ' 

So at 0(1/ Nf) 

, , 1 ^ 2 n a n e n _ / 1 A , 

where r/ is the corresponding critical exponent. Clearly at leading order in 1/Nf, when Nf is 
large, the coefficients {a{\ are accessed. To determine them one needs to compute r\ directly in 
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the 0(1/ Nf) expansion in (f-dimensions. This is the aim of the paper for the operators discussed 
earlier. 



Before reviewing that we make several parenthetical remarks. In assumimg a /3-function of 
the form (2.2) we are restricting ourselves to a particular class of theories which includes QED 
and QCD. If the two loop term of (2.2) had been quadratic and not linear in Nf and likewise the 
three loop term cubic in Nf and so on, then it would not be possible to determine a simple form 
for g c at leading order in 1/Nf. Instead an infinite number of terms of /3(g) would be required. 
This would imply a large Nf expansion would not be possible in that case. This is similar to 
the large A^ expansion of QCD. Then the structure of the /3-function has this nasty form and it 
is not easy to study QCD in a critical 1/N C approach. 

The method to compute critical exponents corresponding to the anomalous dimensions of 
operators in powers of 1/Nf is based on an impressive series of papers, [17-19]. In (l(], [lTj the 
O(N) a model was considered and the technique has been developed for fermion and gauge 



theories more recently, [20, 21]. Essentially one studies the theory precisely at the fixed point 
g c where there are several simplifications. First, at g c there is no mass in the problem so all 
propagators are massless. Second, the structure of the (full) propagators can be written down, 
in the approach to criticality. Therefore in momentum space a fermion and gauge field will have 
the respective propagators ip and A^ of the form, in the limit k 2 — ► oo [16], 



(k 2 y- a ' ^ v ; (k 2 y~p 



(2.5) 



where d = 2/x, b is the covariant gauge parameter and A and B are amplitudes. The dimensions 
a and (3 of the fields (in coordinate space) comprise two pieces. For example, 

a = ii - 1 + \j] (2.6) 

where the first term is the canonical dimension of the fermion determined by ensuring that the 
kinetic term in the action is dimensionless. The second term is the critical exponent correspond- 
ing to the anomalous dimension of ip or the wave function renormalization evaluated at g c . It 
reflects the effect radiative corrections have on the dimension of ift. 

With (2.5) one can analyse any set of Feynman diagrams in the neighbourhood of g c and 
determine the scaling behaviour of the integral. In particular one can examine the 2-point 
Schwinger Dyson equation at criticality to obtain a representation of those equations. It turns 
out that one obtains a set of self-consistent equations which can be solved to determine rj ana- 
lytically as a function of d. Furthermore the approach is systematic in that 0(1/N 2 ) corrections 



can be studied too. In IS] this was extended to n-point Green's function. If, for example, 
one considers the 3-point interaction then the exponent or the vertex anomalous dimension is 
found by computing the (regularized) set of leading order integrals with (2.5). The residue of 
the simple pole of each graph contributes to the anomalous dimension. We will illustrate these 
remarks explicitly in the next section. However, we note that the regularization that is used is 
obtained by replacing (3 of (2.5) by (3 — A. Here A is assumed to be small like the e used in 
dimensional regularization, [ISjl . 

We conclude this section by recalling another feature of critical theory which is important 
in analysing QCD in large Nf. So far the above remarks have been completely general and 
summarize the approach taken in other models. Another common feature is that the theory 
that underlies a fixed point is not necessarily the unique model describing the physics. More 
than one model can be used to determine the (measured) critical exponents. In this case such 
theories are said to be in the same universality class. From a field theoretic point of view one 
can use this to simplify large Nf calculations. For example, the O(N) a model and (ft theory 
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with an 0(N) symmetry are equivalent at the d-dimensional fixed point where the former is 
defined in 2 + e dimensions and the latter in 4 — e, [15]. The critical exponents computed in 
either are the same. For QCD there is also a similar equivalence which has been demonstrated 
by Hasenfratz and Hasenfratz in [[12]. They showed that as Nf —* oo QCD and a non-abelian 
version of the Thirring model are equivalent. The lagrangians of each are, for QCD, 



and 



L 



L 



il 



il 



(G%? 
4e 2 

2A 2 



(2.7) 



(2.1 



1 < I < N r , 1 < a < N? 



for the non-abelian Thirring model, (NATM), where 1 < i < Nf 

l,D fMlJ = d^u + TfjAl, G%, = d^Al - d v A% + f ahc A b ^A c v , f' 1 is the quark field and A® is 
the gluon field. The coupling constants of each lagrangian are e and A and are dimensionless in 
4 and 2 dimensions respectively. The auxiliary spin-1 field of the NATM can be eliminated to 
produce a 4-fermi interaction which is renormalizable in strictly two dimensions. Essentially at 
a fixed point it is the interactions which are important and which suggest (2.7) and (2.8) are 
equivalent. The quadratic terms serve only to define the canonical dimensions of the fields and 
coupling constants as well as establishing various scaling laws for the exponents. Although the 
NATM does not appear to contain the triple and quartic vertices typical of a Yang Mills theory, 
it was demonstrated in j2l| that as Nf — > oo that they correctly emerge from the three and 
four point functions of (2.8) involving a quark loop. Therefore to compute fundamental critical 
exponents in QCD at 0(1/Nf) it is sufficient to use the much simpler theory (2.8) to compute 
these as only one interaction arises. We will point out later where the contributions from 3-gluon 
vertices would arise in the exponent calculation with respect to perturbation theory. Having 
defined the NATM in which we will calculate, we define (3 as 



1 



n 



x 



2.c 



where x is the quark gluon vertex anomalous dimension. Also in using a gauge field in a covariant 
gauge the ghost sector needs to be added to both (2.7) and (2.8). However it turns out that at 
0(1/ Nf) there are no contributions to any exponent we compute and we therefore omit them 
here, ||. 

Finally we state several earlier results which will be required. The /5-function of QCD is 



2 -T(R)N f - ^C 2 (G) 



/3(g) = (d-4)g + 

+ l -C 2 (R)T(R)N f + b -C 2 (G)T(R)Nf 



i7 C 2 (G) 



12 



1 1 7Q 

C 2 (R)T 2 (R)Nf + -j^C 2 (G)T 2 (R)N f 2 



9 

205 



72 



+ 



1 

16 



C 2 (R)T(R)N f 



432 
1415 
864 



C 2 (R)C 2 {G)T(R)Nf 



2857 



Ci(G)T(R)N f + —C}(G) 



g 4 + 0(g 5 ) (2.10) 



where the three loop term was given in [24] and the colour group Casimirs are defined as tv(T a T ) 
= T(R)5 ab , T a T a = C 2 (R) and f acd f hcd = 5 ab C 2 (G). Although in our earlier ansatz we omitted 
a constant term in the one loop coefficient its contribution to g c does not appear until 0(1/N?) 
and so 



T(R)Nf +0 {N^ 



(2.11) 
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Various basic exponents are known to 0(1/Nf) and we note 



[(2» -!)(»- 2) +»b]C 2 (R)r,<? 
?1 " (2 M -1)(m-2)T(*) 

where r, = EZi r]i(e)/Nj, r?° = (2/i - - 2)r(2 / u)/[4r 2 (/i)r(^ + l)r(2 - /x)] and fL6| 

[(2/x-l)(/x-2) +A x&]C 2 (iZK [(2 / u-l) + b(^-l)]C 2 (G)r ? ° 
Xl ■ (2/x-l)( M -2)T(i?) 2(2/x-l)( // -2)T(i2) l " J 

Throughout that paper we will work in an arbitrary covariant gauge. We note that the physical 
operators which occur in the operator product expansion have gauge independent anomalous 
dimensions and by including a non-zero b this will give us a minor check on the corresponding 
exponent calculations. The combination z = A 2 B arises too and 

1 2(2 M -l)(/i-2)T0R) 1 • > 

3 Unpolarized operators. 

We illustrate the large N* technique by computing the critical exponent of the simplest operator 
which arises in the operator product expansion. This is the twist-2 flavour non-singlet operator, 

OjW»-/*n = e-^Sip 1 ^ 1 D» 2 . . . DV"T?jip J -trace terms (3.1) 

where S denotes symmetrization on the Lorentz indices. Although this has already been treated 
in 1/Nf in [14] its value forms part of the flavour singlet calculation detailed later. The full 



critical exponent associated with (3.1) is r/^ , 

r$ = V + r)o (3.2) 

The first piece corresponds in exponent language to the wave function renormalization of the 
constituent fields of (3.1). The second part reflects the renormalization of the operator itself. 
Although each term of (3.2) is gauge dependent the combination is gauge independent. In 
perturbation theory the renormalization is carried out by inserting (3.1) in some Green's function 
and examining its divergence structure. Here we determine the scaling behaviour of the integrals 
where O ns is inserted in a quark 2-point function. The two leading order Feynman diagrams 
are given in fig 1. With the regularization each graph is evaluated with the critical propagators 
(2.5) in cl-dimensions. As in perturbative calculations Q we project the Lorentz indices of the 
operator into a basis using a null vector An, with A 2 = 0. (This is not to be confused with the 
regularizing parameter A which is a scalar object.) They have the general form, omitting the 
external momentum dependence, 

| + Y + 0(A) (3.3) 

where X and Y are functions of d. The integrals are straightforward to compute using standard 
rules for massless integrals. To obtain the leading order large Nf contribution a and (5 are 

replaced by /i and 1 respectively. Following |1| the residue X of each graph contributes to t]q \. 
In this instance we have for the respective graphs 



2liC 2 {R)t,° 



(a* — 2)(2^t — l)T(R) 



1-6 ^ 



(H + n-l)(ji + n-2) 
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and 



( At -2)(2//-l)T(i?)^(/i + /-2) 



E 



i 



(3.4) 



where we have included a factor of 2 in the second to account for the contibution of the mirror 
image and used the value of z, (2.14). Summing the contributions yields, 14], 

2C 2 {R){n - l) 2 ??? r (ra-l)(2// + ra-2) 2/i 



(n) 



(2//-l)( M -2)T(i2) 



+ 



(/i + n- l)(/x + n-2) (ju-l) 



[^(/x + n-1) - ^(/i)] 



(3-5) 

where ^(x) is the logarithmic derivative of the T-function. We recall that this result is in 
agreement with all known perturbative MS results to three loops, ]|, ^]. In concentrating on the 
detail for this operator we will follow the same procedure in the remainder of the paper with 
minimal comment. 

We now turn to the treatment of the flavour singlet twist- 2 operators. Before analysing at 
the fixed point we need to recall several features of their perturbative renormalization. First, 
the operators are, [25, B], 



F 

G 



•n-l 



S^j^D^ 2 . . . D^ip- 7 - trace terms 
\i n - 2 StxG a ^ v D^ . . . D^G 11 Jf n - trace terms 



(3.6) 
(3.7) 



As each operator has the same dimension in four dimensions and quantum numbers they mix 
under renormalization. In other words defining the vector Oi = {Of,Og} then the bare and 
renormalized operators are related by 



bare 



(3.8) 



where 2Ty is a 2 x 2 matrix of renormalization constants. Consequently the associated anomalous 
dimension is a 2 x 2 matrix 'jij(g). It has the following structure, with the Nf dependence 
explicit, 



where, for example, 



91 



jQ9 ^99 



a\g + (a 2 N f + a 3 )g 2 
cyNfg + c 2 N f g 2 



a i 



Cl 



a-2 



2C 2 (R) 



45i(n) 



\(n 2 + n + 2)T(R) 
n{n + l)(n + 2) 

"4 160 



n(n + 1 / 
di = 



big + {b 2 N f + b 3 )g 2 

(diJVJr + d 2 )g + (d 3 N f + d 4 )g 2 



±{n 2 + n + 2)C 2 {R) 



(3.9) 



n(n 2 



-T(R) 



T(R)C 2 (R) 

16[lln 7 + 49n 6 + 5n 5 



32 

() S 1 (n) + -jS 2 (n) 



329n 4 



514n 3 - 350n 2 - 240n - 72] 



9n 3 (n + l) 3 (n + 2) 2 (n- 1) 



32C 2 (R)T(R) 



(n 2 + n + 2) 
(n + l) 2 n(n 2 — 1) 



+ 



Si(n) 



(3.10) 



The remaining entries of (3.9) can be found in, for instance, [l2| and are not important for the 
present situation. To compute Jij(g) the operators are inserted in both quark and gluon 2-point 
functions. Various one loop graphs which occur are illustrated in figs 1 and 2. Clearly from 
(3.9) the Nf dependence is not the same in each term. For example, at g c the Nf dependence of 
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each entry is respectively, 0(l/Nf), 0(1/Nf), O(l) and 0(1). Further, in practical applications 
it is sometimes useful to compute with the operator eigenbasis of (3.9) which simplifies the RGE 
involving Jij(g) and therefore the evolution of the Wilson coefficients of the operator product 
expansion. From (3.9) this leads to the eigenvalues 



a\ + d 2 ± \[M)g 



+ 1 ((a 2 + d 3 )N f + a 3 + d A ± 10=) 9 2 + 0(g 3 ) (3.11) 



where 



Ax 
A 2 



d\Nf 



1 + 



2(g?4 — ai) 



+ 



46i£i 



2JV / [(di(d 3 -o 2 ) + 2ci62)JV / 

+ {da - ax)(d 3 - a 2 ) + <ii(d 4 - a 3 ) + 2(ci& 3 + c 2 &i) 



(3.12) 



Or evaluating at g c the related eigenexponents are at leading order in large Nf 



dxN f g 



ax 



hex 
dx 



g + 



a 2 



b 2Cx\ 2at 



0{Nfg*) 



(3.13) 



Clearly the Nf dependence in each eigenexponent differs. The eigenoperators associated with 
each eigenvalue, A±, are a combination of the original operators. For example, that associ- 
ated with A_ has predominant contributions from the fermionic operator (3.6). Likewise A+ is 
associated primarily with (3.7). 

For the critical point analysis there will be a 2 x 2 matrix of critical exponents analogous to 
(3.9) which are computed by inserting the critical propagators into the diagrams of figs 1 and 2. 
In addition the graphs of fig 3 are also of the same order in 1/Nf. However in determining the 
contribution to X of each of the graphs it turns out that several are trivial due to the imbalance 
of the Nf dependence already mentioned. For instance the leading order term for A + arises 
purely from the tree graph of fig 2. Therefore we take as its entry in rjij = Jij(g c ) as l2~T 



VGG,i 



2c 



(3.14) 



Also rjpQi does not need to be evaluated as its leading order value is given purely by the one 
loop perturbative result. Next the contribution from the final graph of fig 2 is identically zero. 
That is, with (2.5) the graph is A-finite. Therefore the only non-trivial entry to compute is 
ifp-p x- As the non-singlet part has already been determined this reduces to evaluating the two 
loop graphs of fig 3. Each is 6-independent and respectively contribute, for even n, 



m(/x - i)r(n)r(2^° 



(ji -2)(2n-l){fi + n-l)(p + n- 2)T(2fi -l + n) 

x [n(n(n - 2) + 2(p - 2 + n)(2yt - 3) + (2/i - 2 + n)) + 2(n 



2)(p + n- 1) 



and 



Hence, 



(n) 



8/i(/x-i)r(n-i)r(2/i)c 2 (flK 

(/x - 2)(2/i - l)r(2/i - 1 + n)T(R) 



(3.15) 



(2n-l)(fx-2)T(R)N f 



2(jm- l)(n- l)(2fj, + n- 2) 



(M + n-l)0* + n-2) 



+ 4/x[^(/i - 1 + n) - 4>(n)] 



6 



liT(n - l)r(2/j) 



(fi + n - l)(/j + n - 2)r(2/i - 1 + n) 

,2 i „ i o,. q\2 



x [(n 2 + n + 2/i - 2) 2 + 2(p - 2)(n(n - l)(2/i - 3 + 2n) + 2(/i - 1 + n))] 



(3.16) 



We now discuss the structure of rjij. Unlike the perturbative mixing matrix jij(g), rjij 
is triangular. At first sight this would appear to be inconsistent with perturbation theory. 
However, at leading order in 1/Nf the calculation of rjij in fact determines the critical anomalous 
dimensions of the eigenoperators directly. This is not unexpected if one studies the dimensions 
of (3.6) and (3.7) at g c . There clearly the canonical dimensions of each operator is different and 
therefore there is no mixing. The vanishing of certain graphs of fig 2 is merely a reflection of this 
in the large Nf calculation. This indirect relation between the exponents of the eigenoperators 
(3.6) and (3.7) is the reason why we distinguish the perturbative entries of (3.9) by q and g in 
contrast to F and G for the eigenoperators. A further justification of this point of view comes 
from the comparison of the coefficients of the 0(e) and 0(e 2 ) terms in the expansion of (3.16) 
with A_ evaluated to the same order at g c . We have checked that they are in total agreement 
with (3.10) for all n. A further check is that the anomalous dimension must vanish at n = 2. 
Then the original operator corresponds to a conserved physical quantity, the energy momentum 

(2) 

tensor which has zero anomalous dimension. From (3.16) it is easy to check that r/pp , = 0. 
It is worth commenting on this calculation in relation to the NATM and QCD equivalence 



noted earlier, [22]. Clearly A_ and r/pp contain contributions from the insertion of gluonic 
operators in a Green's function. However the graphs we evaluate to obtain rjpp involve only 
(3.6). The resolution of this apparent inconsistency is obtained by studying the integration of 
each quark loop in fig 3 with (2.5) and the e expansion of the individual graphs. Clearly in 
perturbation theory the graphs of fig 1 will contribute to the one loop renormalization whilst 
those of fig 3 will give part of the two loop value of the anomalous dimension. So one would 
expect the large Nf graphs to be 0(e) and 0(e 2 ) respectively. This is not the case. Studying 
(3.15) each graph of fig 3 is 0(e) and from (3.16) their sum is also of this order. The point is 
that after performing the quark loop integral and examining the resulting one loop integral, it 
contains a part which would correspond to the ordinary perturbation theory two loop value as 
well as a piece that corresponds to the final graph of fig 2 which is a one loop integral. In other 
words an effective gluonic operator like (3.7) emerges naturally in the exponent calculation. In 
effect we are confirming in our calculation the equivalence observed in p2| where we recall that 
the three and four point gluon interactions were similarly reproduced by integrating out quark 
loops. 

We conclude this section by giving an indication of the n-dependence of at least the leading 
order 1/Nf coefficients of higher loop terms in the series for 7_ (g). Having established the 
correctness of our expansion at two loops the higher order coefficients are 

6 3 ci 2 10 2 17 2(n 2 + n + 2) 2 [S 2 (n) + S 2 (n)] 

d\ 9 27 27 72 3rr(n + 2)(n + 1) (n — 1) 

2S 1 (n)(lQn 1 + 74n 6 + 181n 5 + 266n 4 + 269n 3 + 230n 2 + 44n - 24) 
9(n + 2) 2 (n + l) 3 (n- l)n 3 
- [100n 10 + 682n 9 + 2079n 8 + 3377n 7 + 3389n 6 + 3545n 5 + 3130n 4 

+ 118n 3 - 940n 2 - 72n + 144]/[27(n + 2) 3 (n + l) 4 n 4 (n - 1)] (3.17) 

and 

6 4 ci 2 „ . , 10 „ . , 2 n . , 2 131 



&4 — — 7— = — SdM — — Sain) — — So(n) — — Si(n) + 

d x 27 y 1 81 v ; 81 y 1 81 w 1296 
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+ C(3) 



2(ra 2 + re + 2) 2 



27 



4(n 2 + 



' 5l(n) 27n(n + l) 9 9n 2 (n + 2)(n + 1) 2 ( 
n + 2) 2 [25 3 H + 35 2 (n)gi(n) + ff 3 (n)] 



27(n + 2)n 2 (n- l)(n + l) 2 

„ r „ . , (16n 7 + 74n 6 + 181n 5 + 266n 4 + 269n 3 + 230n 2 + 44n - 24) 
+ 2 5 2 n + 5f (n p 07 Z^I7 ZTw 7^ " 

- 25i(n)[88n 10 + 608n 9 + 1947n 8 + 3405n 7 + 3670n 6 + 3693n 5 

+ 2973n 4 - 8n 3 - 920n 2 - 48n + 144]/[81(n + 2) 3 (n + l) 4 n 4 (n - 1)] 
+ [68n 13 + 548n 12 + 1861n n + 2474n 10 - 817n 9 - 4143n 8 - 1712n 7 



■ 817n 9 - 4143n 8 - 1712n 7 
2871n 6 - 7702n 5 - 2586n 4 + 3136n 3 + 1952n 2 
288n - 288]/[81(n + 2) 4 (n + l) 5 (n - l)n 5 ] 



(3.18) 



For future reference we list the values of (3.17) calculated for low moments in table 1. A similar 
table was produced for the analogous coefficient in the non-singlet case. It is important to note 
that all the fractions up to n = 8 are in exact agreement with the recent explicit three loop 
singlet results of M, when allowance is made for different coupling constant definitions. 



4 7 5 . 

To apply the large Nf method to polarized operators we need to review the treatment of 7 s 
in perturbation theory and earlier 1/Nf calculations. As is well known one must be careful 
in arbitrary spacetime dimensions when 7 s or the pseudotensor e^ uap are present, The 
simple reason is that both are purely four dimensional objects unlike, say, 7^ and rf" 11 and do 
not generalize in the arbitrary dimensional case. Therefore problems will arise in perturbation 
theory when one uses dimensional regularization. With this regularization calculations are 
performed in d = 4 — 2e dimensions where the infinities are removed before taking the e — * 
limit. There are, however, various ways of incorporating 7 s in such calculations, [26-28]. 
(A review is, for example, [29|.) The original approach of [26] was to split the d-dimensional 



spacetime into physical and unphysical complements. In the former subspace Lorentz indices 
run from 1 to 4 whilst they range over the remaining dimensions in the latter. So, for example, 
the 7-matrices are split into two components 

y = + Y (4.1) 

where the bar, - , denotes the physical four dimensional spacetime and the hat," , the remaining 
{d — 4)-dimensional subspace. Then the Clifford algebra reduces to 

{f\ri = 2tr , {r,r} = o , = &r (4.2) 

The anti-commutativity of 7 5 is not preserved in the full spacetime. Instead the following 
relations are used 

{7 M ,7 5 } = , [f\7 5 ] = (4.3) 

It is known that these definitions give a consistent method for treating 7 s , [28]. Traces involving 
an odd number of 7 s 's are performed via, in our conventions, 

tr( 7 V 7 V70 = ^ vcp (4-4) 

which acts like a projection into the physical dimensions. Further if (4.4) occurs in a loop 
integral where the 7-matrices are contracted with loop momenta the integral is performed first 
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and then the Lorentz index contractions carried out, with the caveat that external momenta are 
physical, = 0. 

For high order perturbative calculations this splitting of the algebra is not always practical, 
p0| . It would be easier if a <i-dimensional calculation could be performed. Such an approach 
has been introduced in [30, and carried out successfully for 3- loop calculations in IpCfl . The 



first step there is to replace 7 5 by 



7 5 = ^W7 M 7V7 P (4.5) 



and remove the e-tensor from the renormalization procedure. The 7-matrices of (4.5) are treated 
as d-dimensional in the calculation before projecting to the physical dimension. If two such e- 
tensors are present then they can be replaced by a sum of products of 77-tensors which is treated 
as d-dimensional. One performs the renormalization in a minimal way as usual to determine the 
renormalization constants. To complete the calculation, in relation to the MS scheme, one must 
introduce a finite renormalization constant Z5 in addition to the first renormalization constant 
in order to restore the Ward identity, J3l[] . 

For the treatment of 7 s in the 1/iVf expansion we recall the simple example of the flavour 



non-singlet axial current. In |33 the method outlined above was followed to correctly determine 
the anomalous dimension of 0ns = ■07^7 5 '0 at 0(1/Nj). First, if one wishes to find the critical 
exponent associated with the non-singlet vector current Ons = tp'y^ip then it is inserted in a 
2-point function and the residue with respect to A is determined. The only relevant graph at 
0(1/Nf) is the first graph of fig 1. If we insert the more general non-singlet operator i/jTip then 
the contribution to the critical exponent from the graph is 

[-f-fTw* " 2/i(l - 6)IK 



2(2 M -l)0u-2)T(i?) 



(4.6) 



where the square brackets are understood to mean the coefficient of the matrix T after all 7- 
matrix manipulations have been performed for an explicit form of T. Therefore for T = 7^, (4.6) 
gives 

[(2/x- !)(//- 2) + buffi 



(2/i-l)(/i-2) 

and so with (2.12) and (2.13) 



(4.7) 



= ( 4 - 8 ) 



consistent with the Ward identity in exponent language, J20|| . For O^s one performs the 7-algebra 
of (4.6) using (4.2), to give 

[(2//-9)(/x -2) + buffi 
(2n-l)(n-2)T(R) [ ■ ' 

Thus 

= (2„ - l)T(R) (410) 

where ~ denotes that the object still has to be augmented by the finite renormalization. As 
discussed in |33f | this does not preserve four dimensional chiral symmetry and is not consistent 
with the Ward identity. To proceed correctly we need to include a finite renormalization constant. 
In this was computed to be 

C 2 (R)e f f ln[l-4T(R)N f a s /(3e)} } /l\ 
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where L is the Laurent operator which removes non-singular terms from the expansion of the 
braces and B(x,y) is the Euler /3-function. The constant Zf s is defined from the requirement 
that, @, 

Zf 1l m $ 0£ 1>) = 7 5 Tl m ® V) (4.12) 

where 'R-j^g denotes the i?-operator or renormalization procedure. In other words the anti- 
commutativity of 7 s is restored by this condition. Using the information in this finite renormal- 
ization together with (4.10) the correct MS anomalous dimension does emerge to all orders in 
the coupling at 0(1/Nf). 

There are several disadvantages, however, with the form of (4.11). First, it is not as compact 
as the 0(1/Nf) exponents that have been produced in earlier work, |16|| . Second by examining 
(4.11) the result can be simplified since the construction of Zf s is in effect equivalent to the 
difference of the exponents (4.7) and (4.9) at 0(1/NA. In other words the contribution from 
the finite renormalization to the final MS exponent is equal to 

8 "° (4.13) 



(2^-l)T(fl) 

Thus the sum of (4.11) and (4.13) correctly gives in MS 

7? m5 = (4.14) 

^ns 



Another difficulty with this procedure is that there is a quicker derivation based on features of the 
fixed point approach. In perturbation theory the regularization used is dimensional in contrast 
to the critical point method. There the spacetime dimension is fixed and the regularization is 
analytic as it is the gluon dimension which is adjusted. The upshot is that, at least for non- 
singlet currects, one can use the anti-commutativity of 7 s in d-dimensions. Therefore with T 
= 7^7 5 in (4.6) anti-commuting 7 s twice immediately gives the same contribution as V = 7^. 
Hence the MS result (4.14) follows directly. We have checked this procedure explicitly for other 
non-singlet operators such as and Sipj^j^ 1 D^ 2 . . . D^tp by calculating the analogous finite 

renormalization constant from a condition similar to (4.12) and observing that the result agrees 
with the direct anti-commuting 7 s calculation. So, for example, the unpolarized and polarized 
non-singlet twist- 2 operators have the same anomalous dimensions, (3.5). In other words we 
have justified the use of an anti-commuting 7 s in non-singlet sectors of calculations. Although 
much of the content of this section may appear straightforward, there is an important lesson in 
the result (4.11) from |^] for singlet operators. Then closed quark loops with an odd number 
of 7 5 matices will occur which means quantities like (4.11) will need to be computed. As we 
have demonstrated that this is equivalent to the difference in the anomalous dimensions of the 
operators of the renormalization condition (4.12) defining the finite renormalization constant, 
flavour singlet operators can be handled in an efficient way. We will come back to this point in 
a later section. 



5 Polarized singlet operators. 

We now extend the unpolarized singlet calculation of section 3 to the polarized case as it is 
important to compare with recent perturbative calculations [9-11]. The twist- 2 operators are 

I, 

Of 1 = i^Stp^j^D^ 2 ...D^ii- trace terms (5.1) 
0P, 01 = \i n - 2 Se^ ia ^trG M D^ ...D^G^ n a - trace terms (5.2) 
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Several features of the computation of the critical exponents will parallel section 3 such as the 
triangularity of the mixing matrix and the Nf dependence of 7? 01 (g) ■ The essential difference is 
the effect 7 5 has in the two two loop graphs of fig 3 which we focus on here. The contribution 
from the graphs of fig 1 is the same as (3.15). 

With (4.4) the second graph of fig 3 is A-finite and gives no contribution to r/pp r For the 
other graph one can compute the quark loop in ci-dimensions before carrying out the second 
loop integral, also in arbitrary dimensions. The projection to four dimensions is made at the 
end. Adding all pieces we have, 



pol 
^FF,i 



2C 2 (R)v? 



(n-l)(2/i + n-l)(/i-l) 



(2/i-l)(/i-2)T(i2) 
+ 2 M (/z-i)[Y>( M -l + n)-^(//)] 



(ji + n- l)(/i + n- 2) 

/j,(2fi + n- 5)(n + 2)T(n)T(2fi) 



(5.3) 



2(ji + n - + n - 2)T(2fi + n - 2) 



As in section 3, due to the Nf dependence we have 



(5.4) 



We have checked that the e-expansion of (5.3) agrees with the anomalous dimension of the pre- 
dominantly fermionic eigenoperator of the mixing matrix at two loops, [9-11]. For completeness 
we note in the notation of (3.9), 



pol 



2C 2 (R) 



45i (n) 



n(n + 1) 



pol 



ol 



8(n-l)T(R) 



jpol 



pol 



6? Gl 



n(n + 1) 
T(R)C 2 (R) 



- 3 T(R) 



^Si(n) + f 5 2 (n) + 



A(n + 2)C 2 (R) 
n(n+ 1) 

32[10n 4 + 17n 3 + 10n 2 + 21n + 9] 



9n 3 (n + 1) 



32(n + 2)C 2 (R)T(R) 
3n(n + 1) 



Si(n) 



8 (n + 2) 
3 + (n + 1) 



(5.5) 



This agreement, moreover, justifies our treatment of 7 5 at the fixed point to be an anti- 
commuting object whose appearance in closed loops is treated with (4.4). Finally we deduce 



a3 



b 3 ci 



pol 



2 10 2 17 2{n + 2){n-l)[S 2 (n) + Sl{n) 
g5 3 (n) - -S 2 (n) - -Mn) + 3n 2 {n + 1 y 



+ 



and 



a>4 — 



-i pol 



2Si(n)(13ra 3 - 6n 2 + 2n + 3)(n + 2) 
9(n + l) 3 n 3 

[61ra 6 + 83ra 5 + 27ra 4 + 217n 3 + 68ra 2 - 36n - 18] 
27(n + l) 4 n 4 



9 If) 2 2 131 

^5 4 (n) - ^s(n) - -5 2( n) - -5 l( n) + _ 



(5.6) 



+ C(3) 



4 2 
27 5l(n) "27n(n + l) 



1 2(n + 2)(n-l) 
9 ~ 



9n 2 (n + l) 2 

4(n + 2)(n - l)[25 3 (n) + 35 2 (n)5i(n) + Sf(n)} 
27n 2 (n + l) 2 
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2(13n 3 - 6n 2 + 2n + 3)(n + 2)[5 2 (w) + Sf(n)] 
+ 27n 3 (n + l) 3 ( ' 

2S"i(n)(49n 5 - 29n 4 + 95n 3 + 41n 2 - 15n - 9)(n + 2) 
81(n + l) 4 n 4 

(19n 8 - 46n 7 + 3ra 6 + 195n 5 - 392n 4 - 317n 3 - 32ra 2 + 54n + 18) 
+ 81(ra+ l) 5 n 5 

In the second column of table 1 we have evaluated (5.6) for low moments as a check for future 
three loop calculations. Although we have given exact fractions for the coefficients at three 
loops, the numerical values of the polarized and unpolarized entries do not differ significantly 
as n increases. 



6 Singlet axial current. 

Having considered a variety of fermionic operators which are both flavour non-singlet and singlet 
we turn to the remaining current. The renormalization of the singlet axial current 0g 5 = 
^7^7 5 ^ is somewhat special. Unlike the singlet vector current the conservation of Og 5 is spoiled 
at the quantum level by the chiral anomaly, [34-36]. Consequently under renormalization the 
composite operator can develop a non-zero anomalous dimension. By contrast the conservation 
of the vector current ensures it has a zero anomalous dimension at all orders in the coupling 
constant. Before attacking the problem of computing the 0(1/Nf) exponent for 0§ 5 in the MS 
scheme, it is worthwhile reviewing the perturbative approach |32| and in particular [[}(]] . (Other 
related contributions to the renormalization of the axial anomaly are [37-39].) In the three loop 
analysis, p0| , two renormalization constants are determined in a manner described earlier for 
other currents. One is the renormalization constant which removes the infinities in the usual way 
but using the standard 7-algebra and the definition of 7 s , (4.5). This renormalization does not 
preserve the axial anomaly, in operator form, in four dimensions. To remedy this a second finite 
renormalization constant Z^ nom is required. The relevant constraint in the present instance is 
determined by ensuring that the operator form of the anomaly, [34-36] , 



is preserved, leading to, [|g, 29], 



z l nom n m ( A W A ) = ^fp n m ( A ^ ap G%G a ap A) 



(6.1) 



(6.2) 



The large Nf calculation follows this two stage approach. In other words the exponent 
corresponding to the MS anomalous dimension of is given by 

Vs = V + V5,s + vl n (6.3) 

It is straightforward to compute the first graph of fig 3 in d-dimensions with the rules given 
previously. With (4.9) 
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{2^ - 1) (/i - 1) 



(6.4) 



where the 6-dependence has cancelled. We have used a split 7-algebra here to be consistent with 
the treatment of closed fermion loops in determining 775 n . There the first graphs of fig 1 and 3 
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will occur as subgraphs. We have checked that the e-expansion of (6.4) agrees with the three 
loop result for the same quantity in [[3C| . 

To compute rj^ n we use the result of section 4. There with a split 7-algebra the finite 
renormalization exponent was determined from the difference in the anomalous dimensions of 
the operators arising in the defining relation. In that case the restoration of the Ward identity 
was simple in that the result obtained was equivalent to using a fully anti-commuting 7 s initially 
and the operators themselves were similar in nature. For 77^ n the exponents of d^O^ and G 
_ e [iuapQa^Qa^^ w \ 1 [ c \ 1 are total derivatives must be determined separately in d-dimensions. In 

detailing that calculation we focus on d^O^ first. 

We insert d^O^ 5 into a gluon 2-point function as illustrated in fig 4. For the moment we 
take the momentum flow to be p into the left gluon leg and (p — q) out through the other. This 
leaves a net flow of q through the operator insertion which is needed since a non-zero momentum 
must contract with 7^7 5 in momentum space. To simplify the calculation of each integral we 
differentiate with respect to q^ and contract with ex^e^P 8 where A and ip are the Lorentz indices 
of the gluon legs. Then q is set to zero, 



This procedure ensures that part of the integrals 
contributing to the renormalization of the operator is projected out. We have given the 0(1/Nj) 
diagrams in figs 5 and 6. The former is the one loop anomaly and with the critical propagators 
it is A-finite. However, as the remaining graphs represent the higher order corrections the value 
of the first graph of fig 5 must be factored off each to leave a formal sum of terms 



6T(R)N f (2fi - l)(/x - 2) 



1 + ^(f + 



(6.5) 



The overall factor is the d-dimensional value of the anomaly which is non-zero in four dimensions 
In (6.5) we have included Z\ from the amplitudes of the quark fields which explains the origin 
of the factor (p, — 2). The residue X is the value of the 0(1/ Nf) part of the dimension of d^O : 
we require. 



/u5 

s 



With the momentum flow as indicated we have computed the value of each graph of fig 6. 
No graphs have been included where the vertex with the external gluon is dressed. These graphs 
together with the vertex counterterm do not contribute to X as they are A-finite in sum. With 
the critical propagators only the first two graphs are non-zero and give 



X 



C 2 (R)y? 

T(R) 



[(2/i-9)(/x-2) + M 
(2 M -l)(M-2) 



+ 



(/*-!) 



(6.6) 



The remaining graphs are each A-finite and we note the colour factors of the last two graphs 
are each C2(G). Recalling the field content of d^Oi? 5 we have 



mo, 1 



T(R) 



(2 M - 1) (ft- 1) 



(6.7) 



The treatment of G is parallel to that just outlined. With the same projection of momenta 
the tree graph of fig 5 gives the normalization value of (— 6) analogous to that of (6.5). The 
relevant graphs are given in fig 7 and we list their respective contributions to X as 

C 2 (R)y? _ [2C 2 (R) - C 2 (G)]r ] ° C 2 (G)[V - 6/x + 1 + %? 
T(R) ' T(R) ' 2(2//- l)(n-2)T(R) 

C 2 (G)[8^ 2 -13M + 4- K^-bM (f R) 

2(2 M -l)(M-2)T(i?) 1 ' j 
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Useful in carrying out this calculation was the symbolic manipulation programme Form, [^Q] . 
The value of the three loop graph accounted for the most tedious part of the calculation. How- 
ever, we made use in part of results of integrals which arose in the computation of the QCD 
/^-function, EJ. This was achieved by computing the dimension of the composite operator 



(G® u ) 2 associated with the coupling constant in a gluon 2-point function. We have included 
a non-zero b to observe its cancellation as a minor calculational check. Although the graphs 
involved in computing the dimension of G in fig 7 are similar in topology to those for dO$ 5 in 
fig 6, the values obtained are somewhat different. For example, the last graphs of each figure 
are similar once the loop integral with the singlet current insertion is performed which leaves a 
Feynman integral with an effective G insertion. The difference in the values arises due to the 
critical propagators used and the fact that this loop integral changes the dimension of the gluon 
lines contracted with it and therefore the nature of the remaining loop integrations. One check 
on this is that the leading terms in the e expansion of each graph ought to agree. It is easy to 
observe that the first two values of (6.8) give the same leading coefficient as the second term of 
(6.6). Likewise the remaining two terms of (6.8) are 0(e). 

With the field content dimension (2.12) and (2.13), we find 

VG,i = - -^y- (6.9) 

It is reassuring to note the cancellation of the terms involving C%(G) again as the overall MS 
renormalization of C?g 5 at 0(1/Nf) is expected to be proportional to C2{R) only. 

With (6.7) and (6.9) the finite renormalization exponent is 



fin _ _ C 2 (R) V ° 

15,1 T (jq 

Therefore 



+ 



3(/i - 2) 



(2/z-l) Ou-1) 



(6.10) 



^ = - („ - l)T(R) (6 - U) 

where the cancellation of the terms proportional to 8/(2//— 1) reflects the non-singlet calculation 
of section 4. A final check on this relatively simple result is that it correctly reproduces the large 
Nf leading order two and three loop MS coefficients of [[}(], p2j1 . This agreement, moreover, again 
strengthens the validity of our treatment of 7 5 . Consequently we deduce, in the notation of (2.1) 
and our coupling constant conventions, 

4G 2 (fi) [9C(3)-7]C 2 (fi) 
«4 = - , «5 = Tl (6.12) 

7 Discussion. 

We conclude our study by remarking on possible future calculations in this area. The natural task 
to be performed next will be the 0(1/ Nf) corrections to the non-singlet twist- 2 operators. Such 
a calculation would mimic the determination of the mass operator dimension but would require 
the quark dimension 772 first. Only the abelian values are available for both these quantities, 



[21]. On another front the corrections to (3.14) and (5.4) are needed. This would parallel the 
calculation of the QCD /^-function in 1/Nf, fUjl . Both these results for the gluonic operators 
would give important insight into the n-dependence of the higher order anomalous dimensions 
and the ^-behaviour of the DGLAP splitting functions. 
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From a more mathematical physics point of view such analyses may become important for 
studying the operator content of strictly four dimensional gauge theories which have (infrared) 
fixed points, [ 42] ) 43]. Evaluating the perturbative anomalous dimension of the composite opera- 
tor at these points would be necessary to gain information on the (conformal) field content of the 
underlying theory in the perturbatively accessible region. Moreover the existence of fixed points 
such as that of Banks and Zaks in QCD, [42], for a range of Nf values have been the subject 
of recent interest in supersymmetric theories with various gauge groups and matter content, 
[43]. Therefore any information that can be determined from traditional field theory methods 
and which sum perturbation theory beyond present low orders such as 1/Nf, could be used 
to compare estimates of, for example, critical exponents deduced from exact non-perturbative 
arguments. 
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A DGLAP splitting functions. 

In this appendix we discuss the relation of our results to the DGLAP splitting functions, P(x, g). 
The x-dependence of the higher order contributions to these functions is currently of interest in 
relation to low x physics. We recall that the anomalous dimensions of the twist-2 operators are 
related to the DGLAP splitting functions via a Mellin transform with the proviso that x, which 
is the variable conjugate to the moment n, is restricted to the unit interval. First, we note that 
our convention is 

dxx n ~ l P{x,g) = \^ n \g) (Al) 

Of course (A.l) can be evaluated at the fixed point and a critical exponent which sums the 
leading order in 1/Nf of the splitting function can be deduced. So, for example, the twist-2 
non-singlet anomalous dimension (3.5) gives the result, 



P ,(xa) - 2 ^)(/^-lW 
*ns,l\x,g c ) 



(2 M -l)(M-2)T(i2) 

- ^-l)x^~ 2 {l- x) + 



^-^ + 1) 5(l-x) 



(p-1) 2 



(/i-l) 



(A.2) 



where we use the standard notation in the final term to ensure sensible behaviour in the x — > 
1 limit. We have checked that the e-expansion of (A.2) correctly reproduces the 0(1/Nf) part 
of the two loop non-singlet splitting function of || . In light of this we determine the three loop 
structure as 

ns _ 17 (l + x 2 )ln 2 x (lbc 2 - 12x + ll)lnx 2(7 - 6a?) 2 

«3 - 72 ^ ^ " 18(1 -x) 27(1 - x) rT~ + 27(1 - x)+ ( 3) 

The treatment of the unpolarized and polarized singlet cases are similar but we detail only 
the latter case as the exponent has a simpler n-dependence. Moreover recent articles have 



examined the polarized splitting functions in the small x limit, [46]. Performing the inverse 
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Mellin transform the splitting function becomes 



p ll 0>5c) 



(2 M -1)(M-2)T(12) 



(/x 2 - 4/i + 1) 



(5(1 -x) - (/i-l)x^ 2 (l-x) 



+ 



2x^~ 2 2x^" 2 
(/"-I) ~ (/i-l)(l-x)+ 
(2/x-l)(/i-3)(^-4)x^ 2 



//(2^-l)x^" 2 (l-x) 



(A.4) 



[xBi_»(2m-2,1- /i) - Bi_ a .(2^-2 J 2-M)] 



(m-i) 

where -B^ (p, g) is the incomplete /3-function which has the integral representation 



B x (p,q) 



(A.5) 



We can once again perform the e-expansion of (A.4) and attempt to compare with the explicit 
two loop results of |u], [ll], [J. We find from (A.4) 



«i 



02 



&1C1 



6 2 ci 



pol 
pol 



6(1 + x) lnx 



- 15x + 8x 2 ) 



(A.6) 



1 c/h \ , H N , 9 „„, H x , /- \ 2(6 — 5x 2 )lnx 

— (5(1 -x) + 1 + x ln 2 x + 10(1 -x lnl-x + 

12 3(1 — x) 



4(l + x)[Li 2 (x)-Li 2 (l)] + 



2(13 - 21x + 13x 2 
9(1 -x) 



(A.7) 



where the dilogarithm function Li 2 (x) enters, [45]. Its appearance, however, would seem to 
suggest that (A.7) does not relate to information in the two loop splitting matrix since Li 2 (x) 
is absent at leading order in 1/Nt there, []l0| , 11]. (It does occur, for example, at next to leading 
order in 1/Nf in all entries bar 7".) The resolution of this rests in the triangularity property 
of the critical point mixing matrix. In order to correctly compare with the explicit perturbative 
matrix it has to be mapped to a similar structure. This is achieved by an invertible 2x2 matrix 
R. In other words, in matrix language, 



P s tri (x, 5c ) = R / dxx"- 1 P s P ert (x, 5c ) )R- X 



(Ai 



As we are only interested in the Oil/NA FF component this reduces to comparing (A.6) and 
(A.7) with the sum of the products of the appropriate elements of (A. 8) in x-space which is a 
combination of entries similar to the left side of (A.6). Also since R is n-dependent we need 
to express this product of n-dependent functions as a Mellin transform of a single x-dependent 
function for the comparison. Useful in this respect is the convolution formula for the Mellin 
transform. As splitting functions are defined to be zero outside the unit interval this takes the 
following form in this instance 



A4[/i(x)]A4[/ 2 (x)] = M 
where the Mellin transform is defined to be 

M[f{x)\ = 



du 



fl 



dxx n ~ l f{x) 



(AS 



(A.10) 
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So with (A. 9) and, at leading order, 



pol 



„P°1 



pol 



C 2 (R) 



1 + x 



(1-x). 



- 2(2x - l)T(R) , d\ 
C 2 (R)T(R) 



pol 



3 
2 

2T(R) 



5(1 - x 

5(1 -x) 



b Po\ 



(2-x)C 2 (i?) 



1 ; (ll-12x + llx 2 ) (l + x 2 ) i 

0(1 — X) + ; : + — ; r In X 

12 v ; 9(1 -x)+ 2(1 -x) 



— (1 — x) + (1 — 3x)lnx + (l + x)ln 2 x 
AC 2 (R)T(R) 



\(x + 4) + J(2-x)ln(l-x) 



(A.ll) 



it is straightforward to verify that (A. 6) and (A. 7) emerge. 

Finally, having established the relation of (A. 4) with perturbation theory it is a simple 
exercise to produce 



A3 



pol 



¥-5(1 -x) - -(1 + x)lnxln 2 (l - x) - — (1 - x) ln 2 (l - x) 
T2 3 3 



— lnxln(l — x) 



9 



ln^x 



(6 - 5x 2 ) 2 

— ; — in x 

9(1 - x) 



4 10 
+ -(l + x)lnx[Li 2 (x)-Li 2 (l)] - — (1 - x) ln(l - x) 
o y 

2(8-9x)(2 + 3x) 2 . in ,. T . . . ... 

+ 27(1 -X) ^ + g( 5 + lla; )[ Ll 2^)- Ll 2(l)] 

+ ^(l + x)[Li 3 (x) -Li 3 (l) -lnxLi 2 (x)] - -Li 2 (l-x) 

8/ n \ tt • / 1 \ W1 w . m (55 - 108x + 55x 2 ) 

+ - 1 + x Li 3 1-x -In 1-x Li 2 1-x + ^ ^ r 

3 27(1 — x)+ 

(A.12) 

In relation to the work of [M we deduce from (A. 6), (A. 7) and (A.12) that the leading small x 
behaviour of these perturbative coefficients are, respectively, 



61nx 



In 2 x 



l ^ 3 
— in x 

9 



(A.13) 



In respect of the large Nf and small x limits, [46] concludes that these do not commute. Although 
this may seem to be a disappointing result it is worth recalling that the primary motivation of 
this paper is the provision of information on higher order coefficients of operator dimensions 
which can be compared with explicit perturbative calculations. 
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n 


^3 - hci/di 


[a 3 - & 3 ci/rfi] po1 


2 





7_ 

1UO 


4 


121259 


3496333 


720000 


19440000 


6 


3166907 


26205857 


13891500 


111132000 


8 


1328467729 


466437737839 




5038848000 


1728324864000 


10 


304337312935261 
1054350180576000 


309708615382541 
1054350180576000 


12 


842357166098254633 
2737572318857376000 


853943993349679513 
2737572318857376000 


14 


42512567719680559 
131614053791220000 


343857572061363287 
1052912430329760000 


16 


755896148277147625515451 

99^1 971 ^^70^4409^407^000 


762758545067630140156811 

99^1 971 ^^70^4409^407^000 
zzjiz i ±uou / yoi^uzo^y ( uuuu 


18 


1121815282809553973842772849 
3235896767484248927963904000 


1130341459796219185897273169 
3235896767484248927963904000 


20 


78640886458671664340562623 
220772683941420492161280000 


26382991363083553371777301 
73590894647140164053760000 


22 


4248342909129791924572989157741 


4271977516708367936099843954621 


11650178316263341224273468364800 


11650178316263341224273468364800 



Table 1. 0(l/Nf) coefficients for unpolarized and polarized fermionic twist-2 singlet 
operators at three loops as a function of moment n. 
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Fig. 1. Leading order graphs for r]^' . 
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Fig. 2. One loop graphs for singlet operators. 
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Fig. 3. Additional graphs for singlet operators. 
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Fig. 4. Operator insertions in gluon 2-point functions. 




Fig. 5. Leading order graphs for d^Og and e lLl " Jp G a iW G a a insertions. 
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Fig. 6. Graphs for d M 0s insertion. 
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Fig. 7. Graphs for e» uap G a G% p insertion. 
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